Effective Theory for Trapped Few-Fermion Systems 
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We apply the general principles of effective field theories to the construction of effective interac- 
tions suitable for few- and many-body calculations in a no-core shell model framework. We calculate 
the spectrum of systems with three and four two-component fermions in a harmonic trap. In the 
unitary limit, we find that three-particle results are within 10% of known semi-analytical values even 
in small model spaces. The method is very general, and can be readily extended to other regimes, 
more particles, different species (e.g., protons and neutrons in nuclear physics), or more-component 
fermions (as well as bosons) . As an illustration, we present calculations of the lowest-energy three- 
fermion states away from the unitary limit and find a possible inversion of parity in the ground state 
in the limit of trap size large compared to the scattering length. Furthermore, we investigate the 
lowest positive-parity states for four fermions, although we are limited by the dimensions we can 
currently handle in this case. 

PACS numbers: 03.75.Ss, 34.20.Cf, 21.60.Cs 
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I. INTRODUCTION 

The properties of strongly interacting Fermi gases have 
been the object of great interest in recent years. Feshbach 
resonances allow the tuning of the interaction between 
trapped particles so that one can study the evolution 
from a dilute Fermi gas to a Bose-Einstein condensate. 
In three-dimensional optical lattices one can reach the 
low-tunneling regime where each site is an essentially iso- 
lated harmonic trap occupied by few fermions [l], [2] ■ This 
opens up a new window into the study of few-body sys- 
tems with two-body scattering lengths 02 that are large 
compared to the range rg of the interaction, \a2 \ ^ tq. 

Large- 1 02! systems are of particular theoretical interest 
when particle momenta Q are small compared to l/rp, 
because then they exhibit universal behavior, that is, the 
system properties depend essentially only on 02 (and for 
bosons or multi-component fermions, also on a three- 
body parameter), but not on the details of the inter- 
action. (For a review, see Ref. Q.) The presence of 
a harmonic trap introduces another parameter, the trap 
frequency uj or, equivalently, the trap length b = l/y/JEJ, 
where ^ is the reduced mass of two particles. As long as 
6 ^ ro, the trapped system should still exhibit universal 
behavior, which for 6< \a2\ could be significantly differ- 
ent from that of the untrapped system. In the unitary 
limit 1 02 1 — > 00, the untrapped two-body system has a 
bound state at zero energy, and a collection of two-state 
fermions is characterized purely by the parameter that 
sets the size of the system. 

While large- 1 02 1 systems have been popular in atomic 
physics mainly in the last decade, they have been inves- 
tigated in nuclear physics since its beginning. The two- 



nucleon (NN) system has two S'-wave channels where 
1 02 1 ^ tq: the scattering lengths are about 5 fm and —20 
fm in the ^Si and ^So channels, respectively, compared 
with a range of about 2 fm. Here we combine two meth- 
ods previously developed to deal with untrapped particles 
— effective field theory (EFT) and no-core shell model 
(NCSM) — in order to provide solutions of few-fermion 
systems in a harmonic trap. 

Since we are interested in the long-range dynamics, we 
can approximate the complicated short-range physics as 
a series of contact interactions that are delta functions 
with an increasing number of derivatives. This can be 
formulated as a non-relativistic EFT where observables 
are expressed in an expansion in powers of Qro (assum- 
ing, for simplicity, that the size of all scattering param- 
eters except for the scattering length is set by ro). The 
EFT in the untrapped two-body sector [1] reproduces 
the effective-range expansion and is equivalent to a pseu- 
dopotential, but can be extended to more-body systems 
[^[^. It can be shown that in leading order the Hamil- 
tonian for two-component fermions consists of a single, 
non-derivative, two-body contact interaction between the 
different components. For systems of identical bosons or 
more-component fermions, a non-derivative, three-body 
contact interaction is also present at this order. The 
EFT for short-range forces and its generalization for the 
exchange of light quanta are reviewed in Ref. 7] . 

The NCSM is a powerful many-body technique for 
solving the Schrodinger equation for A strongly interact- 
ing particles, where the many-body basis states are con- 
structed using harmonic-oscillator (HO) wave functions. 
In nuclear physics, NCSM is used to describe proper- 
ties of light nuclei without adjustable parameters 
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Starting from interactions that describe the NN scat- 
tering phase shifts and selected few-nucleon properties 
with high accuracy, a unitary transformation is used to 
construct effective interactions in truncated spaces of a 
finite number of energy excitations on top of a minimum 
configuration. For the two-body problem, the maximum 
number of excitations is equivalent to the number Nmax 
of shells included. After a truncation where higher-body 
terms in the effective interaction are neglected — the so- 
called cluster approximation — numerical diagonalization 
allows a good description of nuclear energy spectra and 
other properties. Unfortunately, there is no a priori jus- 
tification for the cluster approximation in the context of 
phenomenological interactions. 

We have recently proposed [lOi a combination of 
NCSM and EFT in which the cluster approximation is 
seen as a consequence of the EFT power counting. In- 
stead of performing a unitary transformation on phe- 
nomenological interactions, we simply solve the EFT 
within the truncated space. We determine its parame- 
ters from some binding energies and then predict other 
bound-state properties by extrapolation to the limits 
uj ^ and N^ax oo. We have successfully applied the 
method in leading order to systems with A < 6. How- 
ever, application beyond leading order becomes cumber- 
some because of the increased number of parameters that 
would have to be adjusted to properties of light nuclei. 
This motivates us to devise a more flexible approach, 
which will allow us to determine two-body parameters 
from two-body data. 

In this paper, we present the new method and apply 
it to the case of two-component fermions in a harmonic 
trap. The two [llj and three [12,, T3] -body systems with 
a pseudopotential in a trap have already been studied. 
We provide an alternative, explicit solution for these sys- 
tems in a NCSM basis where u is kept fixed but Nmax 
is made large. This allows us to test the accuracy of 
our approach in a non-trivial system. In addition, we 
present a first solution of the four-body system. The pa- 
per is organized as follows. We introduce our approach 
to the renormalization of the many-body problem in Sec. 
im and present our results in Sec. IIIII Conclusions and 
perspectives for future applications arc discussed in Sec. 

IIYI 



II. THE MANY-BODY PROBLEM AND 
RENORMALIZATION OF THE INTERACTION 

Here we consider a system of A identical two- 
component fermions of mass m, where the two-body scat- 
tering length 02 is large compared to the range rg of the 
interaction, \a2\ ^ ro- This system is assumed to be 
trapped in a harmonic potential of frequency w, whose 
length parameter is also sufficiently large, b ^ ro. We al- 
low various values of the ratio 6/02 ranging from the uni- 
tary limit 6/ 1 02 1 to the untrapped case &/|a2| — 00. 
(Note that we use units such that h = c = 1.) 



Under these conditions, details of the inter-particle po- 
tential are irrelevant, and in leading order can be replaced 
by a non-derivative two-body contact interaction between 
different components. Without loss of generality we can 
refer to the two components as spin-1/2 states, in which 
case the contact interaction acts only in the ^Sq channel. 
We denote its interaction strength by Cq. We start with 
the Hamiltonian for trapped particles. 
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where n and pi represent the position and momentum, 
respectively, of particle i, and [i < j] denotes a pair of 
particles with opposite spin. This Hamiltonian can be 
written as a sum of relative and center-of-mass (CM) 
motion, so that the energy of the system is the sum of 
internal (E) and CM (Ecm) terms: 
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The CM motion is that of a simple HO, so in the following 
we focus on the relative motion. 

The internal Hamiltonian is not well defined as it 
stands, since the delta function is singular. Renormal- 
ization is necessary: after a truncation of the Hilbert 
space with a cutoff in energy or momentum, Cq is taken 
to depend on the cutoff in such a way that finite, cutoff- 
independent results are obtained for observables. The 
many-body problem cannot in general be solved analyt- 
ically, but renormalization can be incorporated easily in 
numerical calculations, which are formulated from the 
outset within a finite space. 

These issues can be made explicit in the two-body sys- 
tem, where the relative motion of two particles of oppo- 
site spin is described by 
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in terms of the relative coordinate f* and the reduced mass 
fi = m/2. 

In the conventional NCSM approach, one chooses the 
model space by truncating in the number of shells in- 
cluded in the calculation. More precisely, one consid- 
ers only HO states with the principal quantum number 
N = 2n+l smaller than a given value N„iax- The effective 
Hamiltonian is then constructed via a unitary transfor- 
mation U {H2^^ — UH2U^), designed so that the lowest 
D eigenvalues of H2 are exactly reproduced by H2^-^ {D is 
the dimension of the model space). Of course, the trans- 
formed Hamiltonian no longer has the form of the initial 
interaction, and, in particular, additional non-local terms 
are induced even if one starts with a contact interaction. 

The alternative approach proposed in Ref . [l^ is based 
on principles of EFT, and it constructs the interaction 
between particles making use of power counting. Thus, 
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in each model space one preserves the form of the inter- 
action, as dictated by the power counting, and one deter- 
mines the strength of each interaction strength as a func- 
tion of Nmax and oj so that some observables are exactly 
reproduced. For example, one can fix Co — Co{Nmax,'-^) 
such that the two-body ground-state energy is fixed. 

In this paper we propose an intermediate approach. 
Thus, as in Ref. m we consider only the terms dic- 
tated by power counting, but instead of adjusting their 
strength to reproduce observables in the few-body sys- 
tem, we reproduce the w-dependent energies of two-body 
states, as in the conventional NCSM approach. The num- 
ber of states whose energies are to be reproduced is fixed 
by the number of coupling constants to be adjusted, and 
not by the dimension of the model space, as in the con- 
ventional approach involving a unitary transformation. 
We note that, according to the general principles of EFT, 
even if initially we start with only a contact interaction, 
each truncation is going to induce additional correction 
terms. Effectively, the truncation of the space induces 
effective range, shape and other parameters, which can 
be adjusted to the appropriate values as more derivatives 
of delta functions are included. 

To renormalize the interaction, we consider the relative 
Hamiltonian ^ and solve the corresponding Schrodinger 
equation. 
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in a finite model space defined by < 2n < N^ax- Since 
the delta function acts only on S waves, we consider here 
only those waves. The solution ip{'r) can be expanded in 
the complete set of HO wave functions 
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with (j)n{f) the 5-statc solutions of the unperturbed HO, 
If 2n\ 



4)n{r) = 
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in terms of generalized Laguerre polynomials. One can 
easily verify that the eigenvalues of Eq. (|3|) are given by 
the consistency condition 
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At this point, the coupling constant C^iN^riax-,^) is 
undetermined. We can fix the value of C^iNmaxT'^) at 
each Njnax by requiring that we reproduce a given mea- 
sured level E{uj) of two trapped fermions 0. However, 
the experimental values are very close to the the theo- 
retical prediction of a pseudopotential, given [ll[ by the 
transcendental equation 



r{3/A-E/2uj) _ b 
T{l/4:~ E/2uj) ~ 2^' 



(8) 



where a2 is the two-body scattering length in the absence 
of the trap. (For discussions about the limits of appli- 
cability of Eq. (HI), see Ref. Q). The solutions of Eq. 
([5]) come in levels, and we can use any of the low-energy 
states to fix the value of Co(-/Vmaa;, w). Here we use the 
most natural choice, the ground state. 

Once we adjust the strength of the two-body term for 
each Nmax value (and, in general, to), we can perform 
few- and many-body calculations. The many-body model 
spaces for A > 2 are chosen so that the two-body trun- 
cation is included consistently in relative coordinates 0] . 
Note that if one wants negative-parity states, one has 
to truncate the many-body space to odd Nmax- In this 
case, we use a two-body interaction that is adjusted to 
the largest even number below N^ax, that is, N„iax — 1- 

Results for the energies of few-fermion systems are in 
general complicated functions of N^ax and lu. In the 
unitary regime the dependence on uj gets simpler. When 
6/ 1 02 1 0, the two-body spectrum of Eq. ([8]) is given 
by En/oj = 1/2 + 2n, where n = 0, 1, . . .. Taking into 
account the b dependence of Eq. we see from Eq. 
^ that the cu dependence factorizes: Co{Nmax,<-^) = 
{2Trb/ fi)joiNmax), where ^o{Nmax) is dimensionless and 
uj independent. Indeed, in this case both /z and b dis- 
appear from Eq. ^ upon a rescaling r ^ bp (with p a 
dimensionless variable) , and the energy can only be pro- 
portional to u!. More generally, rescaling all coordinates 
in Eq. ^ will ensure that ui exactly factors out in all the 
many-body eigenenergies. This was to be expected, given 
that for 1 02! 00 the only energy parameter is w, which 
also sets the size of the system. It has been argued 
that the proportionality constant between E and u is re- 
lated to the short-distance scaling exponent 7 of the wave 
function. The three-fermion spectrum has been solved 
semi-analytically in Ref. [T^ , while the ground-state en- 
ergies for A = 2 ^ 22 have recently been calculated in 
the case of a short-range potential . 

In the opposite limit 6/|a2| — > 00, the trap is removed. 
Eq- & gives expected results. There is a single two- 
body bound state of energy Eq = — l/2/za|, if the two- 
body interaction is sufficiently strong, 6/02 — > 00. All 
other states have energies En/u! = —1/2 + 2n, where 
n — 1,2, . . ., and correspond to scattering states. When 
b/a2 — > — 00, the interaction is weak and the spectra of 
all few-body systems are expected to approach those of 
non-interacting fermions in a HO. This can be seen from 
Eq. ([7]): as 6/02 ^ —00 the two-body ground-state en- 
ergy approaches 3w/2, leading to a divergence from the 
n = contribution and to Co{Nmax, — > 0. In this limit 
all few-body energies should be set by occupation of HO 
levels. For example, all other two-body levels are given 
by the poles in the right-hand side of Eq. ([7]). Whether 
few-body bound states exist when the interaction is suf- 
ficiently attractive will be discussed below. 

In any case, at the end of the calculation we extrap- 
olate to the limit Nmax 00. As we show in the next 
section, this extrapolation is relatively smooth. Our ap- 
proach can be straightforwardly extended to higher or- 
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ders. At the two-body level, we add to Eq. ((T|) in- 
teractions proportional to derivatives of the delta func- 
tion, starting with an S'-wave range correction of relative 
OiQro) [J. A P-wave interaction (the most important 
between like components) appears at 0{{Qro)^). For 
spin- 1/2 fcrmions contact three- and more-body forces 
act in P and higher waves, so these forces are of relative 
0{{Qro)^) or higher Many-body properties are thus 
to a high degree determined by two-body forces alone. 
The two-body parameters can be determined from an 
extension of Eq. ([8| to non-negligible-range interactions 



III. RESULTS 

In this section, we apply the method introduced in the 
previous section to the description of systems of three 
and four identical fermions in a harmonic trap. We con- 
sider the unitary regime (|a2| — > oo) as well as a general 
non-vanishing 6/02 value, for both positive and negative 
scattering lengths. 

In the three-body calculation, we have employed two 
numerical methods. In the first stage, we have used the 
three-body code in relative coordinates employed in Ref. 
[To| . Thus, to handle three identical fermions of spin 
1/2, we allow interaction only in the ^Sq channel, calcu- 
lating the isospin T = 3/2 solution. In a second stage, we 
have developed a new program that can handle fermions 
of arbitrary spin, dropping the isospin from the possible 
quantum numbers. Both codes produce the same results 
for fermions of spin 1/2, but the latter can be extended 
to larger Nmax- We should point out that in order to 
correctly antisymmetrize the three-body system, we in- 
clude all possible I states in each model space. There- 
fore, although we show below only results for the levels 
that are affected by the contact interaction, we also ob- 
tain energy levels that are left unchanged. (For example, 
states with j'^ — can be obtained from addition of 
two I = 1 relative angular momenta and total spin 1/2 
or 3/2; however, in such channels the contact interaction 
is not present since it appears only for ^ = 0.) 

Figure [T] presents the main results of our investigation 
of the unitary regime. We show the dependence of the 
lowest-energy levels (labeled by their spin and parity, j'^) 
in units of uo as function of the ultraviolet cutoff Nj^ax- 
We have verified numerically that, indeed, to factorizes, 
as we have argued in the previous section. The running of 
the low-lying energies is fairly smooth when one increases 
the ultraviolet cutoff, and allows an extrapolation to the 
limit N„iax 00. However, we find that the running 
depends on the state considered. In order to qualify this 
statement, we assume a running of the form 
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{Nmax + 3/2)" ' 



(9) 



with -Boo, Ec, and a fitting parameters. In a previous 
publication [lO| , we assumed a running for the energy of 
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FIG. 1: (Color online) Dependence on the ultraviolet cutoff 
Nmax of selected low-lying energies (in units of the HO fre- 
quency) of three harmonically trapped fermions of spin 1/2 
in the unitary regime (6/02 = 0). We present the lowest three 
states of j'^ = ^ (continuous curve), j'^ = i (dotted curve), 



and j" 



3 + 



(dashed curve). 



the many-body system of the form Eq -\- Ei/A, where 
A = \/2^{Nmax + 3/2)0; was defined as the ultravio- 
let cutoff. (Because here the solutions depend trivially 
on uj, we do not have to include it in the definition of 
the cutoff.) This was motivated by the running in the 
continuum two-body system. Here we carry out a more 
detailed investigation of the running of the three-body 
solution with the ultraviolet cutoff. Allowing a to vary 
in the fit, we find values between 0.5 and 2.0. We find 
that, in general, the values of E'oo are not very sensitive 
to the value of a, but can depend on the features of the 
running and the values of Nmax included in the fit. For 
example, while in Fig. [1] the first j'^ = ^ excited state 
appears to be smooth, a closer examination reveals some 
structure, as shown in Fig. [2l In this figure we exhibit 
also three possibilities for fitting, where we include all 
points, only points with Nmax > 15, or only points with 
Nmax > 21. Since Eq. ([9]) is valid for large Nmax, our 
best guess for a is obtained in the case when we con- 
sider only Nmax > 21. In this case we find a ~ 1.474, 
and the extrapolated value for the energy of the state is 
Eoo/u) ~ 4.83. 

In Tables Ul IHII we present our estimates for Eoo for the 
states shown in Fig. [TJ In our approach, there are two 
sources of errors: on one hand, we neglect higher-order 
terms in the expansion of the interaction (terms involv- 
ing derivatives of contact interactions) and, on the other 
hand, we perform a fit of the form of Eq. ^ to obtain the 
large- iVmaa; limit. The variation with Nmax allows us to 
estimate errors associated with missing terms; thus, for 
low-lying states, even fairly small values of Nmax produce 
results within about 10% of the extrapolated values. On 
the other hand, we find that, in general, the extrapola- 
tion errors coming from fits to large values of Nmax are 
very small. For the levels shown in Tables Ul lIIII these 
extrapolation errors are beyond the last digit shown, ex- 
cept for second level in Table HI where the error is 2 in 
the last digit. 
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FIG. 2: (Color online) Running of the first | excited state 
(circles), and three fits with the functional form in Eq. (|9]) 
taking into account all points (dotted line), only points with 
Nmax > 15 (dashed line) and only points with Nmax > 21 
(dot-dashed line). 



TABLE I: Comparison between the results of the present ap- 
proach (Eoo/i^) and of the semi-analytical formula from Ref. 
[a (Eq. m)Jorr = ^~. 
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q s 


Eq. ^ Eao/uj 





1 


1.77 


2.77 2.76 





1 


1 1.77 


4.77 4.71 
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1 


4.36 


5.36 5.39 



In the limit Nmax —^oo,we should recover the semi- 
analytical solutions obtained in Ref. [l^ under the as- 
sumption of a two-body pseudopotential: 



E 



1 + S + 2q, 



(10) 



where q = 0,1,..., and s > 1.77 is the (real) solution 
of a transcendental equation. The agreement between 
our numerical method and the semi-analytical values, 
also shown in Tables Hl lIIIl is remarkable and provides 
a confirmation of the reliability of our method. As a 
side result, we note that the short-range scaling exponent 
we obtain for three particles is 7 ~ —0.24. (While 
the virial theorem [l^ is not satisfied necessarily in each 
model space, we expect it to be satisfied in the large- 
Nmax hmit.) 

The successful test in the unitary regime encourages 
us to extend the application of our method for arbitrary 
values of the 6/02 ratio. Experimentally, one can vary 
the scattering length of trapped particles by means of a 



TABLE II: Same as in Table U but for f = 
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TABLE III: Same as in Table H] but for j'" = f^. 
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FIG. 3: (Color online) Ground-state energies (in units of the 
HO frequency) of the three-fermion system in a harmonic trap 
as function of the cutoff Nmax for 6/02 = (circles), b/a2 — 
— 1 (squares) and b/a2 = 1 (diamonds). 



magnetic field, and, in principle, obtain a large range of 
ratios 6/02- Thus, for illustration, we have considered a 
range of 6/02 values, and fixed the coupling constant in 
each model space so that the ground-state energy given 
by Eq. ^ is always exactly reproduced. 

In Fig. [3] we show the running of the lowest j"^ = 
^ (ground state) energy of the trapped three-fermion 
system for b/a2 = ±1 in comparison with the unitary 
limit. We obtain different shifts of the energy level with 
respect to the unitary solution depending on the sign of 
02. If 02 > 0, the level is pushed downward in energy 
with respect to the unitary value, while if 02 < the 
level is pushed upward. This behavior is a reflection of 
the same type of behavior in the two-body system: while 
in the unitary limit the two-body ground state is Eq /uj = 
1/2, the corresponding energy for 6/02 = 1 decreases to 
Eq/u! ~ —0.34, and for 6/02 = —1 increases to Eq/uj ~ 
0.9. 

In Fig. m we consider both = i and = 
states, and plot the lowest-energy level in each case as a 
function of 6/02. We show only the extrapolated values 
Eoo obtained by means of Eq. ([9]), taking into account 
only the points with Nmax > 21 to reduce fitting errors. 
For 6/a2^1.5, the ground state is the same as in the 
unitary limit. However, we observe a tendency for an in- 
version of the parity of the ground state for 6/02 > 1.5, al- 
though the positive- and negative-parity levels are nearly 
degenerate. Note that in the region with 02 > we might 
have larger errors from the fitting procedure, because the 
interaction strength increases with increasing values of 
the ratio b/a^. 

As discussed in the previous section, for weak two- 
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FIG. 4; (Color online) Lowest (extrapolated) energy levels 
for j'" — i (circles) and j'^ — i (squares) states of the 
three-fermion system in a harmonic trap as function of the 
the ratio 6/a2. Around 6/02 ~ 1.5, we observe an inversion 
of the parity of the system's ground state. The dashed line is 
a fit of the positive-parity points by a quadratic form. 



body attraction the spectrum should approach that of 
three non-interacting fermions in a HO trap — which is 
itself wide. In the limit 6/02 —00, the lowest state is a 
configuration with two particles in the first S state, and 
the third in the first P level. Figure 2] indeed suggests the 
asymptotic value Auj (Sw) for the energy of the negative 
(positive) -parity state. 

On the other hand, for strong two-body attraction the 
relative effect of the contact interaction increases. In 
order to understand the possible parity inversion, we fit 
the a2 > results for E/uj in Fig. |4] to a quadratic 
function of 6/02 (dashed line). In the untrapped 6/02 ^ 
00 limit we find that the lowest positive- and negative- 
parity levels are very close, and their energy is 
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(11) 



Thus, within errors, the system of three fermions is at the 
threshold for the scattering of one particle on the bound 
state of the other two (dimer): the system is not bound, 
or very weakly bound, and can be viewed as composed 
of a dimer and an additional particle. The lack of a deep 
three-body bound state follows in fact from a very naive 
argument. The existence of a shallow two-body bound 
state is a consequence of a balance between the renormal- 
ized delta-function attraction and the kinetic repulsion. 
The addition of a third spin- 1/2 fcrmion roughly doubles 
both, so we do not expect a collapse of the three-body 
system — as one does when the number of pair interac- 
tions grows faster than the number of particles, e.g. for 
identical bosons and multi-component fermions. 

The parity inversion of the ground state is therefore 
plausible. In the h/a2 00 limit, the positive-parity 
state likely represents the untrapped 5-wave particle- 
dimer scattering state that one expects to dominate suffi- 
ciently close to threshold. This is consistent with results 
[l3| for near-threshold particle-dimer scattering. 

We have also solved the four-fermion problem in the 
trap. For now, we are limited to smaller model spaces 
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FIG. 5: (Color online) Same as Fig. |3] but for the f = 0+ 
states of the four-fermion system. 



than for three fermions because we use the redstick 
Slater determinant code [18j that preserves only the third 
component of the angular-momentum projection, mixing 
states with good angular momentum. The dimension 
of the many-body basis thus increases significantly. Be- 
cause we truncate the many-body configurations allow- 
ing only a maximum number of excitations on top of the 
minimum solution, we can eliminate spurious CM contri- 
butions when we compute the spectrum of the intrinsic 
motion. (The calculations in relative coordinates and in a 
properly truncated Slater determinant basis are, in fact, 
equivalent.) Moreover, because the Hamiltonian is rota- 
tionally invariant, we obtain eigenstates with good an- 
gular momentum (for four particles, the lowest positive- 
parity state has j = 0), even though most of the indi- 
vidual Slater determinants that compose the many-body 
basis do not have good j. 

Figure [S] exhibits the running of the lowest positive- 
parity levels of four harmonically trapped spin-1/2 
fermions, with the same parameters as in Fig. [31 6/02 = 
0, ±1. Because of the present limitation to just a few val- 
ues of Nmax for the four-body system, our fitting proce- 
dure is likely to be subject to larger extrapolation errors, 
although Fig. [5] shows small variation with the ultravi- 
olet cutoff. Thus, in the unitary limit, if in Eq. ^ we 
use a — 1/2, we estimate Eoo/lu ^ 3.58, while if we use 
a = 1, we obtain Eoo/^ 3.85. (The small number of 
points does not allow us to obtain a reliably from the fit.) 
Both values are in line with a previous calculation of the 
four-body system in the unitary limit 16]. From these 
results we can estimate the short-range scaling coefficient 
7 [iSl for four particles between —0.65 and —0.92. 

At non- vanishing 6/02 we observe the same type of dis- 
placement in the lowest positive-parity level as for two- 
and three-fermion systems. However, our technical lim- 
itations in the four-body system do not allow us yet to 
make a more detailed investigation of a large range of 
ratios 6/02- As expected, the system goes to the non- 
interacting case in the limit 6/02 — > —00, that is, the 
ground-state energy is 1301/2. As for three particles, 
we can extrapolate to the other limit 6/02 — > 00 with a 
quadratic form in 6/02 > 0. We find the binding energy 
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in this limit to be about 40% below the threshold of scat- 
tering of two dimers. However, because our results are 
subject to large errors, we leave a more thorough analysis 
of four-body spectra in different regimes for future work 
using relative coordinates, when it should be possible to 
achieve larger ultraviolet cutoffs. 

IV. CONCLUSIONS AND OUTLOOK 

We have presented a first application of the NCSM to 
the description of three- and four-fermion systems in a 
harmonic trap. Based on the general principles of EFT, 
we have proposed a method to renormalize the two-body 
interaction to be used for few- and many-body calcula- 
tions in finite model spaces. Tests of our results against 
semi-analytical calculations in the unitary regime suggest 
good accuracy in estimation of the energy levels of the 
three-fermion system. Moreover, we have extended the 
application of our method to arbitrary 6/02 ratios. 

There are, of course, ways to improve our results. In 
principle, even if a pseudopotential is a good approxi- 
mation, in finite model spaces we can add a range cor- 
rection to the contact term. This should improve the 
running of observables with the ultraviolet cutoff. Al- 
though in many cases the running is quite small, we could 
improve the calculations for 02 > 0, where the interac- 
tion is stronger and requires larger model spaces to con- 
verge. (In the four-nucleon system we can follow the evo- 
lution of the contribution of different configurations. We 
find that the components with the lowest-energy config- 
uration dominate for 02 < 0, while other configurations 
reshuffle small contributions in each model space. How- 
ever, this situation changes, with a larger contribution 
of the higher shells, for a stronger interaction such as for 
02 > 0. This is not surprising since a stronger interaction 
becomes more important across shells.) Adding higher- 



order terms should decrease the magnitude of Ec in Eq. 
([^ , thereby allowing for a more precise determination of 
observables in few- and many-body systems. 

In principle, we could extend our method even fur- 
ther. First, we can perform similar calculations of other 
observables (e.g., r.m.s. radii), with more particles, as 
well as for fermions of different spin. Especially for the 
latter it will be interesting to investigate the ability of 
the method to predict Efimov states in a finite model 
space. In the case of untrapped multi-state fermions, the 
EFT power counting requires a contact three-body force 
to prevent the collapse of the three-body system Q. 

Finally, the same type of approach can be considered in 
nuclear physics, where addition of a CM harmonic term 
to the intrinsic Hamiltonian produces effectively the same 
type of two-body relative Hamiltonian as in Eq. ^ [l^ . 
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